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ABSTRACT 


The Generalized Rotor Performance (GRP) program is a 
computer program designed for calculating forward flight 
performance of a helicopter rotor system at a specific 
flight condition. It can be used to evaluate either an 
articulated or a hingeless single rotor system in forward 
flight or in a wind-tunnel test. The program was originally 
designed by the Sikorsky Aircraft Company and purchased by 
the United States Navy. 

The goals of this thesis were (1) to reinvestigate the 
theory and logic used in the program, (2) to add selected 
desirable features to the program, (3) to produce a much 
needed Users' Manual, and (4) to run an analysis comparing 
the program's calculated results against manufacturer's 
data. These goals were accomplished and the results of the 
analysis indicated that the program produces highly accurate 
results within the normal cruise range of a modern 


helicepter. 
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I. INTRODUCTION 

The Generalized Rotor Performance (GRP) program is a 
computer program designed for calculating forward flight 
performance of a helicopter rotor system at a specific 
flight condition. It can be used to evaluate either an 
articulated or a hingeless single rotor system in forward 
flight or in a wind-tunnel test. The GRP series of programs 
were originally designed by the Sikorsky Aircraft 
Corporation with the cooperation of the United Aircraft 
Research Laboratories. Work began on this series of 
programs in 1955 and has continued to date. The United 
States Navy purchased this program in 1964 and has used it 
in the Naval Air Systems Command as a helicopter performance 
computer program. 

The goals of this thesis were (1) to reinvestigate the 
theory and logic used in the program, (2) to ensure that the 
Navy's version of this program performed the calculations 
according to the design theory, (3) to produce a much needed 
Users' Manual, (4) to add certain desirable features to and 
correct errors found in the program, (5) to run an analysis 
comparing the progran's calculated results against 
manufacturer's data and (6) to document areas that may need 
further attention. The work was done with the help and 
cooperation of both faculty at the Naval Postgraduate School 
and personnel in the Naval Air System Command. 

The output available from the program includes rotor 
shaft horsepower, rotor profile horsepower, wain fotor 
torque, lift, drag, thrust, H force, rolling and pitching 
moments and other forces and moments calculated in tha 
Shaft, control and relative wind ceference axis systems. 
Also available are the azimuthal historias of the bladets 
flapping angles, rates and accelerations, plis azimuthal and 


radial histories of angles of attack, CL, CD, Mach number, 
sweep angles and air load distribution. In addition, the 
program will calculate a Fourier coefficient series for 
flapping angle, radial station air leads and azimuthal 2 
force distribution. 

The program can be divided figuratively into three main 
routines. The first routine datermines a steady state 
flapping solution. The second routine determines the forces 
and moments generated by this flapping solution. The third 
routine compares the calculated results of lift, drag and 
other options to those desired by the user. If the 
calculated vaives area not within a predetermined tolerance, 
the program will, by use of a first order Taylor series, 
generate new values to reenter into the blades flapping 
routine, 

The program uses a combined blade element/momentun strip 
theory. The blade'ts radius is divided into a maximum of 15 
segments. The blade is advanced around the azimuth ina 
prescribed number of degrees. At each particular azimuthal 
position, velocities are computed in the perpendicular (UP), 
tangential (UT) and radial (UR) directfons. From these 
velocities and a known pitch angl@ distribution, the lecal 
angles of attack are calculated. From these angles of 
attack, the forces at each radial blade segment are 
determined. Once the forces ar@ known, the local moments 
can be found. From the summation of these moments at a 
particular azimuthal position, the blade flapping 
acceleration is calculated. This acceleration, combined 
With the calculated flapping angle and rate, is used to 
advance the blade to the next azimuthal position. A steady 
state flapping condition is assumed tc exist when the blade 
flapping angle and rate at the PSI equal zero and 360 degree 
azimuthal position are within a prescribed tolerance of each 


other. 


This method of calculation accounts for retreating blade 
stall, the reverse flow region and compressibility effects. 
The program can accommodate a full range of geometric and 
design variables including flapping hinge offset, elastic 
flapping hirge restraint, first and second harmonic cyclic 
inputs and spanwise variations in blade twist, local mass 
densities, chord and tip sweep. The program uses no snmall 
angle assumptions and has no restrictions on tif speed, 
forward velocity or advance ratio. The rotor system can be 
oriented in any direction in space and can be given any 
rotor shaft or aircraft roll, pitch or yaw angular 
velocities. The program will perform calculaticns in 
sttaight and level flight or a uniform induced velocity may 
be added to simulate climbs and descents. The GRP will 
accept up to five different airfoil data decks plus one 
blade spar characteristics data deck. Lift and drag 
information is entered into the program in the form of CL 
and CD versus Angle of Attack Tables for up to 15 different 
Mach numbers. The user has a choice of six methods of 
soluticn depending upon the desired restrictions. The user 
has available nine printout options, two of which are 
program debugging options. There are also error printout 
messages that will assist the user having difficulty with 
the progran. 

In the method described above, the flapping motion 
determined is about a flapping hinge offset, and with the 
elimination of all assumptions that the flapping angles are 
small, the rotor control axis (axis of no feathering) and 
the tip path plane axis (axis of no first harmonic flapping) 
ace no longer considered convenient for reference in the 
analysis of rotor blade motion. Therefore, the axis 
selected for use in the analysis is the rotor shaft axis 
system. All forces, moments and angles are referred to the 
shaft axis system with the exception of some of the final 
output forces which are referred to the relative wind and 
control axis systems. 
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With the use of the computer, the GRP program eliminated 
many of the simpliflying assumpticns of the earlier 
Classical theory originated by Wheatley and Bailey in Ref. 1 
and 2. The assumptions of the classical theory did not 
impose serious limitations in low speed flight, but as 
helicopter speeds increase, the inaccuracies inherent in the 
theory seriously limit the usefulness of this method. Tha 
assumptions of the Wheatley-Bailey theory which are not 
Fresent in the GRP program include the following. 

1. The flapping and inflow angles are assuned 

small. 

2. The lift and drag coefficients are approximated 
by a linear and a quadratic variation, 
respectively, with angle of attack. 

3. The effects of Mach number on CL and CD are not 
considered. 

4, Sectional characteristics in the reverse flow 
region are the same as those in the 
conventional flow. 

5. The sectional characteristics of blade twist, 
tip sweep, flapping hinge offset and root cut 
out are ignored. 

While the GRP progran represents a refined approach to 
the classical theory, there are still some basic assumpticns 
which make the GRP subject to error. These are: 

1. Steady state two-dimensional airfoil data are 

used. 

2. Quasi-static blade analysis is used. 

3. The rotational speed about the shaft axis is 
constant. There is no lead cr lag motion in 
the program. 

4. The rotor blade is assumed rigid in bending and 
torsion. 


12. 


5. The rotor,inflow is assumed uniform unless 

varied by the uset. 

6. Spanwise flow is incorporated into the 

calculation of blade angle of attack; however, 
it is assumed that the flow at one segment does 
net affect the flow at any other segment. 

It is felt that ‘the errors induced by the above 
assumptions are relatively small in the normal cruise speed 
region of a modern helicopter. This is the region from just 
below the minimum power airspeed to the maximum ailowable 
cruise speed. The program can not calculate rotor hcver 
power. This is because one of the factors in the 
denominator in the main routine which estimates new flapring 
routine reentry parameters is the advance ratio. Since the 
advance ratio is zero in a hover, the computer would attempt 
to divide here by the nanmber zero. Highly accurate results 
are not available in the airspeed region from hover to just 
below the minimun power airspeed while using the normal 
uniform inflow assumption. This is aregion of highly 
non-uniform flow. Variable inflow may be inputed by the 
user. However, at this time, no information is available on 
how successful this technique is in accurately estimating 
the required rotor system hcrsepower. while the program's 
required technique for induciag harmonic variable flow is 
cumbersome, Bramwell in Chapter Four of Ref. 6 describes 
several methods which could be incorporated into the 
pregram. 

The quasi-static blade analysis assumption is valid 
except in the very high speed region where a Large 
percentage of the retreating blade is in a stalied 
condition. This is an area of aerodynamic hysteresis that 
is influenced by unsteady aerodynamics. While encountering 
lift hysteresis, the amount of lift change above the steady 
state CL is about the same as below. However, the lift 
distribution on the rotor disc will change. Stall is 


delayed and occurs at a slightly later azimuthal statica 
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then predicted by quasi-static analysis. While difficult, 
complicated and computer-time-consuming procedures can be 
taken to reduce these assumptions there is no guarantee that 
the accuracy of the solution will improve. It is felt that 
the present program produces highly accurate results in the 


flight range of interest in a modern helicopter. 


TI. METHOD OF ANALYSIS 


a ee ee er ee 


A. GENERAL COMMENTS 


The GRP uses a combined blade element/momentum strip 
theory in its calcultions. The blade radius is divided into 
a maximum of 15 segments or strips. Figure 1 illustrates a 
typical blade element. In crder to cbhtain a steady state 
flapping solution the program must calculate the time 
history of the flapping motion. This necessitates a 
complete knowledge of the flapping angles, trates and 
accelerations. The time history of this motion must be 
found by solving the highly non-linear equation for the 
summation of moments about the flapping hinge. Figure 2 
illustrates the forces which produce moments on the rotor 
system. The blade has aercdynamic forces acting upwards, 
blade weight acting downwards, centrifugal ferce acting 
outwards, an optional elastic hinge restraint acting to 
return the blade to a zero flapping angle and an inertia 
force resisting any change in blade flapping. The moment 
summation equation about the flapping hinge say be written 


. 


as 


(1) Mange * AY t My + Mee + Hy = Oe 


The moment due to inertia force can be expressed as 16. 
waere Iy is the tlade moment of inertia about the flapping 
hinge. Substituting this into equation 1 and solving for 6, 
the following governing equation is obtained. 


ae 
(2) @ = Mpees + My t Mee + Mane 
ty 


If the angle of attack and local Mach nugber at é6ach 


blade segment were known, the CL and CD could be obtained 


j= E/R 


ba 


Figure 1 Blade Element Diagram 
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Figure 2 Rotor Blade Mcment Diagran 
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from the inputed airfoil infcrmation. The local aerodynamic 
forces and moments can be calculated and the flaprfing 
acceleration determined. Equation 2 can be rewritten as 


ae Ke segments 
(3) Gn = EF ( aMasao, + GMwe + UM: + aMeue: } 
121 Iy 


The flapping angle and rate at the (N + 1) azimuthal 
position could be obtained from the following expressions. 


(4) Baar = Bn *~ en at 
2 
e oe t 
(5) Bas = Bn + @, at + 8. = 


However, since the flapping is periodic in nature and 
has a direct relationship to the azimuthal angle, PSI, the 
values for flapping are solved with respects to PSI, vice 
time. The values of ..1, y and time are related by the 
equation AW =mat. Therefore 


» . < 
(6) g = 38 = 1.08 = 28 
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The governing equation for the flappiny motion now 
becomes 


(8) S) = Mayeo + Mw + Mee + Meng 
ep 
» 


The flapping solution is based on the assumption that 
the angle of attack is known. However, it is not and the 
program must proceed through an iterative process in order 
to determine the inflow ratio, collective pitch and cyclic 
input angles required to generate the desired forces, 
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B. ANGLE OF ATTACK CALCULATIONS 


A very important and basic part of this program is the 
procedure by which the local angles of attack are 
calculated. While the program will calculate angle of 
attack with any angular velocity applied either to the rotor 
system or the helicopter, the development here will describe 
level flight only. The classical approach ignores radial 
flow, UR, and the angle of attack would re obtained as shown 
in Pigure 1. However, as the blade rotates abcut the shaft, 
it will encounter a large variation in radial flow. The 
program attempts to compensate for this radial flow in the 
following manner. ‘Instead of the inflow angle PHI equalling 
the arctangent of UP/UT, it is seat equal to UP/UTUR where 
UTUR is the resultant velocity in the tangential and radial 
direction. This is illustrated in Figure 3. The fitch 
angle is also reduced by the cosine cf the sweep angle. The 
angle of attack is now calculated in the sweep plane. This 
three-dimensional angle of attack is lower than the 
classical two-dimensional angle. 

The program enters the CL, CD tables with this sweep 
plane angle of attack and the sweep plane resultant Mach 
number. The program computes the forces using the velocites 
in the sweep plane, UP and UTUR, and the blade chord 
geometry in the normal plane. Once the forces are computed 
in the sweep plane they are resolved into their respective 
directional forces. 

This three-dimensional angle of attack, due to sweep, 
will delay stall on the rector by reducing the angle of 
attack, This describes what actually occurs on the blade. 
However, it is felt by previous personnel using the progran 
that there is a point where the sweep becomes so large that 
it tends to wash out the lift being produced. The progran 
has been modified to reduce CL by one half for sweep angles 
petween 60 to 72.5 degrees and reduce CL to zero for sweep 
angles between 72.5 and 90 degrees. These high sweep angles 
occur normally only on the inboard tlade segments of the 
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Figure 3 Spanwise Flow Diagram 
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retreating blade in and near the reverse flow region. 

UP, UT, UR and pitch angle for level flight are shown in 
Figures 4 and 5. They are calculated in the following 
Manner. The radial velocity, UR, is calculated as 


(9) UR = (Voosd,) cos¥ cos@ 


The tangential velocity, UT, has two components, The 
first is the local rotational velocity, ac, the second isa 
sinusodial component of the forward flight velocity. The 
general expression for UT is 


(10) UT = Sur + (Voosds) sin 
This expression contains a small angle assumption in the 
terme for blade flapping angle. The program acccunts for 


the fact that the flapping angle does reduce the true radius 
slightly by using the following formula 


(11) UT = (E/R + (cr - E/R)cosB)M + (Vecosas)sin¥ 
The perpendicular velocity, UP, consists of three terms, 


the inflow ratio, a flapping velocity and a small component 
of forward velocity. The inflow ratio is defined as 


12 = Ysinws - 
(12) y= 


The second component is a verical flapping velocity 
which is a function of flapping rate and radius. This is 
computed as 


® 
(13) UP(2) = (r - B/R)B 
The thirc component is due to the fact that there is a 


stall component of the radial flow which acts in the UP 
direction due to blade flap angle. This is equal to 


eo. 


Cb) 
[we —}+——_( R~ E/R) cos 8 —_{ 


Y_- (V cos ds) cosy 
tc) Neosds 


(N cos vs) SINY 


UT = (E/R + (rv - E/R)cos® lait (Veosds)ainY 
Figure 4 UT Diagram 


21. 


1 
ee Seem sak tes atts ye abbtk8SUSa ony avid alata it Cada enlace RB Sse Se se cuecen pall CApatbatalass ad bed ak aN“ Mat al ae ak 


Ur = Roos? + (r - E/RIe + (Veosd,; )cos¥sing 


UR = (Vcosdg)cos¥cos@ 
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(Vcosas) cos’ 
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Figure 5 UP and UR Diagram 


(14) 0P(3) = (Vcosa)cos¥ sin’ 
The total formula for UP is 
(15) UP = A Roose + (r= £/R)6 + (Vcosds) cos¥sin® 


The pitch angle (98) is expressed in equation 16. 4.75 
is the pitch angle at the 75 percant racius position, 6 is 
the twist depending on the relationship of the location of 
the blade segment to the cr = .75R location. The next four 
terms are the first and second harmonics of cyclic pitch and 
(tan@,)@ is the coupled effect of the flapping angle on 
the pitch angle. 


(16) @ = (O@a¢ + 8 = AlScos¥ — B1SsinY¥ - A2Scos2¥ 


~ B2ssin2y - tan(%,) @ ) COS aves 


The local Mach number is calculated as Uya where ais 
the speed of sound and U is given in eguation 17. 


‘ 
(17) 0 = (up* + ur? + opty! 


The angle of attack can now be calculated. Figure 1 
illustrates thatet =o-¢g@ . Initially, the program requires 
estimated values for the inflow ratio, collective pitch at 
the PSI equal zero position, harmonic cyclic inputs, blade 
twist and initial flapping angle and rate. The user can 
either input these values or accept the program's automatic 
default values of <-.02, 5, -1.2, 7.53, OO, O and QO 
respectively. Using these assummed values, the initial 
angle of attack can be determined as fcllovws. 

UP = AMRcosh + (cr - B/R)@ + (Vcos dg) cos#sin® 
reduces to 
uP = ALR 
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UT = (E/R + (rc ~E/R)cos#).2 + (Vcosds)sin¥ 
reduces to 

UT = (E/R + (rc -E/R)) SL 

UR = (V¥cosds) costcos® 


teduces to 
UR = Vooseds 
In the program VooSe, is replaced by the term wOR where 

is the advance ratio in the shaft axis system. Since the 
local pitch distribution and inflow ratio are estimated, the 
local angles of attack can be determined. The next section 
describes the method used for calculating flapping angle and 
rates at the N + 1 azimuthal position. 


Sei atiltatinl sah cota Lite dh tris ash Manat siete ceaziasd haw ee ; 


C. METHOD OF SOLUTION FOR THE FLAPPING EQUATION 


In the preceding section, equation 8 was developed in 
order to calculate the time history of the flapping moticn. 


*% 
(8) @ = Macro * Swit ee * “ene 
Ty no 


The relation involves a complicated second order 
differential equation for establishing the flapping angle as 
a function of PSI. The numberical solution is accomplished 
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by use of a finite difference equaticn and a step-by-step 
procedure, An important characteristic of the solution is 
that it is periodic in nature. Tw: function which 
represents a steady state flapping solution has the property 
Gy) = @ (wx2w). Using this fact, a Fourier harmonic series 
can be written to describe the blade flapying mction. 


(18) @ = AQ - Alcos¥ - Bisin® - A2cos2¥ - E2sin2¥ - 
A3cos3¥ - B3sin3¥ .... 


By assuming that the first harmonic flapping is much larger 
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than the other higher harmonics, the series can pe reduced 
to 

(19) @ = a0 - Alcos¥ - Bisin¥ 

This equation can be aifferentiated with respect to PSI 
to obtain 

(20) @ = Atsin¥ - Bicos¥ 

(24) “B= atcosw - Bisin¥ . “3 

Assuming that the values of @, @ ana @ are known at 
some azimuthal position, the following equations must hold 

(22) 6. = AO - Alcos¥® ~- Bisinw 

(23) e. = Alsiny - Bicos#¥ 

(24) “@, = Alcos¥ + 31sin¥ 

These equations can be solved for the N+ 1 azimuthal 
position by substituting Ya = Y%e*s¥ into the above 
formulas. The flapping angle equation becomes 


(25) @ = a0 - Atcos(Hh +89) = B1sin( Yq +A) 


By using the following two identities equation (25) can 
be rewritten as equation (28). 


(26) cos(Y, +Ay) cos¥, cOSAY - sin®& sinay 


\ 


(27) sin(®& *a%) = sin cosaAy + cos% sina’ 


(28) Gna 7 AO ~ Al(cos¥cosay -sin® sina) 
- Bi(sin%cosay - cos® sinay) 


The terms can he rearranged into equation 29. The same 
*« 
procedure can be used to develop equation 30 for g: 


(29) Gna, = AD ~ cosae(Atcostn + B1sin Ws) 
“ + sinaw(Alsin¥n - Bicos'm) 
(30) Ona = COS av(Aisin a - Bicos Ya) 
+ sinaw(Aicos\ + Bisin sy) 


Substitution into equations 22, 23 and 24 reduces these two 
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expressions to the flapping equations used in the progran. 
The user can either enter the value for AW or accept the 
programs automatic default value of 15 degrees. 
” 4 —% 
BN Ban = GB, Cos he + Bn SIN ON 
a 2% 
(3 2) Bee = n+ i ay « (\- cos AY) Bn 


While this integration scheme is not one of the standard 
methods used, it is very useful in obtaining periodic 
solutions for differential equations similar to the one used 
here. Notice that for small values of SY the trignometric 
expression can be reduced to an ordinary Taylor series. By 
assuming sin4SW equals 4¥Y and cosAW equals one eguation 29 
reduces to 

+** 


% * 
(33) Bay = Bn * S* Pr 


By assuming sin AW equals AY and ccsA equals the first 
2 
two terms of the cosine series 1 - AY equation 32 can be 
reduced to 


4 < _** 
34 By, = 8B. + ay@, + 5 CANNY 6, 


D. FORCES, NOMENTS AND RADIAL INTEGRATION 


The forces acting on a rotor blade may be found by the 
summation of the elementry forces along the span at any 
azimuthal position. The forces considered by the progran 
are the resulting aerodynamic forces only, and are initially 
summed in the shaft reference axis system. The fprogram 
computes forces in three axis systems. They are the (1) 
shaft, (2) control and (3) relative wind axis systems. 

The program radially integrates differently for lift and 
drag calculations. The drag is calculated from the hinge 
offset, E/R, to the rotor tip. The Lift is calculated from 
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the hinge offset to the next to last rotor blade segment. 
The last segment is considered a blade "Tip Loss Factor" 
segment. It is assumed that the tip trailing edge vortices 
cause no lift to be produced in this segment. The normal 
procedure is to define this segment as the last three 
percent of the rotor blade radius. 

The first of the maximum 15 hblade segments is 
considered the spar or cut out segment. This is defined as 
the area between the hinge cffset and the point where the 
airfoil actually begins. [If no spar data are entered, it is 
asssumed *hat this first segment produces no lift and drag 
is obtained by using the CD verse Alpha Tables for the first 
blade airfoil data deck entered. If spar data are entered, 
the first segment and all other segments designated spar 
segments will have the lift and drag characteristics of the 
entered spar data. 


E. AZIMUTHAL INTEGRATION METHOD 


Onve the integration of the rotor forces and moments 
along the blade is complete, an integration around the 
azimuth must be performed in order to obtain the average 
forces and moments. Since the soluticn of the flapfing 
equation is obtained by a step- by-step method, the 
integrands of the integrals cover YW are known only at a 
certain number of equally spaced points around the azimuth. 
For equilibrium flapping the integral is a periodic function 
of  , and for this case integration can be shown tc be 
equivalent to an averaging process. This result makes the 
azimuthal integration simple. The following method is used 
where N is the number of azimuthal positions used in the 
calculations and b is the number of fLiades. 
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(39) 


(40) 


Once the program has calculated 


desired solution 
this solution? 

certain known or assumed values 
ratio (A), collective 
B1S, 
3, is to iterate on 
modification of A and 8.75, where 6.75 is 


resultant are shown in Figure 7. 


MAJOR ITERATION 


a steady state flapping 


and has determined the resulting forces and 
a question remains to be answered. Is this the 
and, if not, what must be done to obtain 

In order to solve the flapping equations 


were used, including inflow 
pitch (8.75) and the 
A2S and B2S). One method, 
the required lift and drag threugh 
the pitch angle 


eyclic pitch 
first described in 


at the 75 percent radius station at the PSI equal zero 

azimuthal position. The modified values are then reentered 

into the flapping routine and the calculation is repeated 

until it converges to within a specified tolerance on the 

required lift and drag. Drag is used here in the sense of 

negative rotor propulsive force. The program procedure is 
‘ outlined below. 

The required rotor lift and propulsive force are 
expressed in terms of the magnitude and direction of the 
resultant force in the longitudinal plane. These are shown 
in Figure 6 where it can be seen that 

(36) a" = og + at 

2 2 in 
(37) Re = (Fe) . (Fxg) ) 
b= i ' = ) 
(38) Ry COs (a") D Rp Sin(a") 
The required lift and propulsive force and their 


In a similar way, 


if 
% 5 
R = (lg ¢ Dg ) 
a" = arctan (De/%e) 
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The differences between tha required and the computed Ry 
and a“ values are defined as 


(41) A Ry = Ryo > RL 


"1 " - " 
(42) Aa ag a 
In order to correct \ and 6.75 to compensate for the 
difference between (Ree o'er) and (Rp at) , the required 
values are expanded in a Taylor series with A and 8.75 as 
variables. The first order equations are: 


BRa YWe 
(43) Ree = Re + Sa AY te Jax AS8.s5 


nals Yen + Ya" 
(44) Ap =z A+ A AOS 
: aan 28as 


Solving the equations for the iteration on NN and @.75 


yields 
aa" pQy — sRe aa! 
(45) AA = 2225 2 Ons 


ea 
Wa va" _ 3s@a dal 


Rs Jo" _ 3@e da" 
2 Ons 3075 ON 


Now that Roos Rg and ane and a" are known, the corrected 
values of A and @.75 can be approximated by: 


(47) Awe nN ~ aa Sasn = Sas ~ Nas 


In order to solve equations 45 and 46, the values of the 


partial derivatives in these equations must be found. The 
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>: Figure 6 Resultant GRP Force Calculation Diagram 


Figure 7 Required GRP Force Diagram 
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procedure used in based on the Wheatley-Bailey method and 
the formulas can be found on pages 186 and 207 of Ref. 4. 
Reference 5 outlines the derivation and a complete 
derivation was performed and verified in conjunction with ; 
this thesis. : 

; 


G. FLOW CHART 


The flow chart in Figure 8 is a klock diagram showing 
the relationship between the various parts of the GkP 
progran. 
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III. GRP USERS" MANUAL 
The GRP currently has the capability to enter over 200 

individual case variables and option selectors. Presently, 
170 are available to the users. This manual explains the 
input and output format of this GRP program. Examples of 
both input and output are given in order to make the use of 
this program easier. This manual is divided into the 
following areas of discussion. 

A. Main Iteration Options 

B. GRP Data Deck Order 

C. Airfoil Data Deck Requirements 

D. Spar Data Deck Reguiresents 

E. Sample Data Input Format 

F. Case Input Listings 

G. Case Input Default Values 

H. Case Input Data 

I. Case Input Format 

J. Case Optional Output Indicators 

K. IBM 360 Execution Control Cards 

L. Sample Program Output 

It is recommended that the user examine the GRP Case 

Input listing carefully, since a few options require certain 
variables to be inputed which are not necessarily lecated in 
the same general area of the input listings. An attempt has 
keen made to list all input variables concerned with the 
different options in the discussion cf each option and input 
variable. 


A. MAIN ITERATION OPTIONS 


The GRP offers the user six different options for 
determining a solution. They are as follows, 


t. Normal Routine 
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In the normal solution to the problem the computer 
will vary ‘the inflow ratio, LAMBDA, and the pitch angle, 
@.75, in order to produce the required lift and drag. The 
cyclic inputs, A1S and B1S, are considered constants. 
Uniform inflow is assumed unless the user induces a 
non-uniform inflow by the use of variables 117 and 118, LAML 
and UVL. The program will calculate the required shaft axis 
angle and position of the control axis from the fixed value 
of B1S. In all of the options, the user has a choice of 
using either a program calculated first estimate of flapping 
angle, velocity and acceleration or an initial value of zero 
for all of the above. Either method usually requires the 
Same approximate amount of computer time. The variable 
PCNV, item 97, determines the method to be used. 


2. Desired Flapping Angles 


This option allows the user tc specify the desired 
longitundinal and iateral flapping angles with respect to 
the rotor shaft axis. The program at each intermedient 
force iteration will vary A1S, B1S, inflow ratio and pitch 
angle. Variables 100 and 112, 113 and 114 control the use 
of this option. 


3. Short Iteration Scheme 
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This option follows the normal routine with one 
major exception. The program will make only one fass 
through the flapping routine on each major iteration. The 
program may or may not arrive at a steady state solution for 
flapping in the first few iterations. It is assumed that 
the first few iterations in the normal routine are only 
rough estimates of the way the variables should be changed 
and that an exact flapping solution is not actually required 
if the user is only interested in transient flapping 
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behavior. In order to use this routine, the user must input 
a negative number for the variable XITLIM, item 73. The 
absolute value of XITLIM will still determine the maximum 
allowable number of times the program will enter the major 
iteration routine, 


This option will follow the Normal method described 
in paragraph one. The Normal method will only converge on 
lift and drag and will not consider the moments produced. 
If the variable TRIM, item 160, were assigned a non-zero 
value, the program would attempt to trin out the pitching 
and rolling moments about the rotor shaft. It is suggested 
that this be done in a two case run. The first case should 
be a normal run, with the desired printout. Then, for the 
second case, input the variable TRIM. This will do two 
things. First, it will provide a ccnverged solution without 
consideration of moments. Secondly, it will allow the full 
number of iterations to be used to reduce the moments. The 
program does this by a short routine varying A1S and B1S. 
During this process, the whele flapping and force iterations 
must be repeated, but it will at least start by using a 
converged solution. The variable PCNV, item 97, should be 
non-zero to allow the flapping solution from the previous 
case to be used as a first estimate in this case. Setting 
the variable SKIPIN, item 91, equal to zero will allcwa 
force and moment summation to be outputed for each major 
iteration. This will allow the user to see exactly how the 
program is proceeding. 


9+ [POP _ Option 


In this option, the program iterates upon the 
required lift but ignores any values inputed for required 
drag. This option can be used to simulate a wind-tunnel 
test. The user must input the shaft angle, item 117, and 


cS 


the pitch angle, item 87. These two inputs will be held 
constant. The program will iterate on the inflow ratio, 
LAMBDA, in order to obtain a solution. The variable TOP, 
item 96, controls the use of this opticn. 


6. ALOPT Option 


This option, like TOP, is a wind tunnel option. It 
also iterates on trequired lift only. However, here the 
shaft angle, item 111, and the inflow ratio, item 88, are 
required inputs and are held as constants. The program will 
iterate on the pitch angle, item 87, in order to determine a 
solution. The variable ALOPT, item 110, controls the use of 
this option, 


B. GRP DATA DECK ORDER 


Data is entered into this program in the following 
order. 
1. Airfoil Lift Coefficient Table 
2. Airfoil Drag Coefficient Table 
(Repeat steps one and two as necessary.) 
>. Case Irput Data 
4. Harmonics of, the Inflow Ratio ** 
5. Spar Lift Coefficient Table at 
6. Spar Drag Coefficient Table ** 
** Optional Data 


C. AIRFOIL BLADE DECK REQUIREMENTS 


The GRP program requires that all CL and CD information 
be entered into the program in tabular form. Tables for up 
to 15 different Mach numbers and five different airfoils can 
be entered. The program currently does not use cr require 
values for the Coefficient of Moment, CM. Since certain 
segments of the retreating hlade are in the reverse flow 


region, angle of attack tables are required to include 
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valued from -180 degrees to +180 degrees. If they are not 
included, an error message will be printed when the program 
can not locate a value for CL and CD at these large positive 
and negative angles of attack. If complete angle of attack 
information is not available for a particular airfoil, the 
user can use the values provided in Section E. It is 
realized that available data on airfoil behavior at large 
angles of attack are very limited, but so is the region on 
the rotor disc where the blade operates at these high 
angles. Since this occurs only immediately arcund and 
within the reverse flow region where dynamic pressure is 
low, little precision is lost in performance calculation by 
using one common representation for most airfoil behavior. 

As a minimum, two values for CL and CD 2t two different 
angles of attack and Mach numbers must be supplied. As a 
maximum, 15 different values of Mach number may be entered, 
Each Mach number may contain up to 48 different values of CL 
and 48 different values of CD and their associated angles of 
attack. 

The first Mach number must be equal to zero. This table 
can be an exact duplicate of the lowest Mach number table 
the user has available. The highest Mach number table 
should be high enough in crder to prevent the program from 
stopping because of a local Mach number higher than that in 
the table. A quick check can be obtained by adding together 
the rotor tip velocity and the forward flight speed. This 
combined velocity, dividea by the local speed of sound, must 
be less than the maximum Mach number entered into the 
program. The program linearly interpolates between Mach 
numbers and angles of attack in order to determine the value 
of CL and CD. The subroutine BLINY does the interpolaticn, 

Several options are available to help reduce the number 
of data points that must be entered. If the airfoil is 
symmetrical, the user only needs to enter values for 
positive angles of attack. The program will assign the 
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appropiate sign to the value of CL and CD according to the 
sign of the angle of attack. This is accomplished by case 
input variable 107, SYM. This option can also be used for a 
cambered airfoil where values of CL and CD at negative 
angles of attack are unknown. 

Values for large angles of attack need not be entered 
for each Mach number by making use of the program's input 
variables 156 and 157, or their automatic default values. 
Values for large positive and negative angles of attack need 
only be entered for the two lowest Mach number tables. The 
lowest Mach number table must be at a Mach number equal to 
zero. If an angle of attack is greater than variable 156 
(HIALFA) or Lower than variable 157 (LOALFA), or the 
programs default values of plus and minus 30 degrees, the 
Mach number is set equal to zero. This ensures that cnly 
the first two Mach numbers have to carry the whcle range of 
angles of attack from ~180 to +180 degrees for a cambered 
airfoil or from zero to +180 degrees for a symmetrical 
airfoil. 

The format of the table input will now tke described. 
The first data card contains the variable WELADE. This 
controls whether or not the user receives an echo printout 
of the Blade Airfoil Data being entered. If the value of 
WBLADE is equal to zero, the user will not receive an echo 
printout of the Blade Data. If WBLADE is a non-zero number, 
the user will receive the echo printcut. The read format 
for WBLADE is F10.0. WBLADE is the only item on the first 
data card. 

The next card also contains only one piece of data, 
NBLADE. NBLADZ is the number of different airfoil data sets 
to be used and appears on this card in an I2 format. This 
program will accept up to five different blade airfoil data 
sets. It will also accept one blade spar data set. If the 
rotor blade being analyzed were composed of thrse different 
types of airfoils, NBLADE weuld equal three. However, if 
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the blade consisted of three sections, of which the first 
and third section were the same, NBLADE would equal two. [+t 
is explained later how the blade segments are assigned their 
respective airfoil type. This arrangement provides the user 
with the ability to vary the make-up of the blade while only 
having to enter into the program once a Darticular set of 
airfoil data. 

The above two variables, WBLADE and NBLADE, are only 
entered once for each complete computer run which uses the 
same set of airfoil data. The following information will be 
entered twice for each type of airfoil used. It will be 
entered first for CL and secondly for CD for each type 
airfoil. The overall format is summarized below. 


Card 1 WBLADE 
Card 2 NBLADE 
Card 3- CL's for airfoil number one 


CD's for airfoil number one 
CL's for airfcil number two 
CD's for airfcil number two 
Repeat as necessary. 

Card number three contains the variable NZ, which is the 
number of Mach numbers for which CL‘'s will be entered for 
the first airfoil. This number must te right justified in 
I2 format. The maximum number of Mach numbers for each CL 
and CD for one airfoil is 15. This is the only number 
entered on this card. 

Card number four begins the actual Mach number, CL 
versus angle of attack tables. The format in this paragraph 
must be repeated for each Mach number. This first card is 
divided inte 11 fields (12, 10F7.0). The first field is a 
two-digit, cright-justified integer in 12 format. {It is 
egual to twice the number of data pairs for this Mach number 
plus two. This tells the computer how many numbers are 
required to be entered for this particular Mach number. A 
data pair consist of one angle of attack and its associated 
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CL or CD. The remaining ten fields are each seven columns 
long in floating point or F7.0 format. These fields begin 
in cclumns 3, 10, 17, 24, 31, 38, 45, 52, 59, and 66. On 
this first card, the field that starts in column three 
contains the number of data pairs at this Mach numher. The 
field that begins in columns 10 is the actual Mach number. 
The remaining eight fields on this card are for the first 
four data pairs starting with the lowest angle of attack and 
increasing towards the highest angle of attack. If a 
symmetrical airfoil option is used, the lowest angle of 
attack is zero. If a non-symmetrical airfoil is used, th: 
lowest value for the first two Mach numbers should he -180 
degrees and for all the remaining Mach numbers the value of 
LOMACH entered or ~30 degrees. 

All the remaining cards for this particular Mach number 
will contain the data pairs. These cards contain ten fields 
each, the first field consisiting of nine cclumns and the 
remaining nine fields consist of seven columns beginning in 
cclumn ten and following the same format as the first card. 
The format is (F9.0, 9F7.0). This card is repeated as often 
as needed. Columns 73-80 are not read and may be used for 
ccmments. 

This procedure is repeated until all the cCL's are 
entered. Once completed the program is ready to enter the 
values for CD. The whole procedure is repeated again, 
Starting with the value for NZ representing the number of 
Mach numbers for which CD's will be entered. If more than 
one airfoil data deck is to be used, the abcve procedure 
will start over again by reading in the value of NZ for the 
Number of Mach numbers to be entered for values of CL for 
tue second airfoil. The number of data pairs for each Mach 
mumber must be between 2 and 48, The number of "twice the 
data pairs plus two" must be between 6 and 98. 

Prior to entering the airfoil data, the user should 


review thea following input variables. 
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1. SYM (107) - Symmetrical and nonsymmetrical 
aicfoil data input ccentrol. 

2. HIALFA (156) - The highest angle of attack for 
which values of CL and CD will be found at all 
Mach numbers. 

3. LOAFLA (157) - The lowest angle of attack for 
which values of CL and CD will be found at all 
Mach numbers. 

4. SPAR (103) <- Number of segments using spar 
airfoil data. 

5. TIPSWP (158) - Amount of tip sweep in degrees. 

6. TPSWST (159) - Blade segment number at which 
the tip sweep begins. 

7. BSPL (120) - Input control variable for spar 
data. 

8. RB(I) (161-175) - Controls the blade segment 
airfoii data assignment. 


GD. SPAR DATA DECK REQUIREMENTS 


The format for spar data are similar to that of the 
airfoil data with the exception that only one set of spar 
data can be entered into the program. Before the program 
will read spar data, input variable number 120, BSPL, must 
have a non-zero value. In addition, variable 103, SPAR, 
must indicate the number of blade seguents which are using 
the spar data. The program automatically assumes that the 
first segment is a spar segment. This is further explained 
in the Case Input section. The spar data are the last tc be 
entered into the program. This is an optional input and is 
uot required. If spar data are not inputed, the progran 
will assume that the one automatic spar section creates no 
lift and has the drag characteristics of the first airfoil 


section entered. 


The format for inputing spar data are as fol. ows. The 
first card contains the variable WRSPAR in F10 0 format. 
A non-zero value of WRSPAR causes an echo printout of Spar 
data to occur. The remaining spar data are handled the 
exact same way as the airfoil section data, starting with 
the variable NZ. There is no input similar to that of the 
airfoit section stating how many different spar data decks 
are being entered since only one is allowed. As before, a 
minimum of two Mach numbers are required to be entered. The 
blade segment printout indicates spar segments by the use of 
a "0" for that segment. 

Input variables associated with SPAR data are as 
follows. 

1. SPAR (103) ~ The number of blade segments using 

Spar data. 

2. BSPL (120) - Input variable which controls the 

input and use of spar data. 


E. SAMPLE BLADE INPOT 


The next several pages illustrates the sample fcrmat for 
a blade which has the following characteristics: (1) an 
echo printout is not required, (2) there are two airfoil 
decks to be read in and (3) the first airfoil deck has nine 
values of Mach numbers for which CL's are to be entered. 
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80-82 


GRP CASE INPUTS 


DESCRIPTION 
Tip speed 
Radius 

Speed of Sound 
Air Density 
No, of Blades 
Forward Speed 


Offset Ratio of 
Flap Hinge (e/2) 


Delta X 

Local Twist 

Local Mass Density 
Local Chord 

Delta PST 


Flap Iteration Limit 


Initial Beta 
Initial Beta * 
Initial Beta ** 


Lift and Drag, 
Iteration Limit 


Required Lift 
Required Drag 
Lift Tolerance 
Drag Tolerance 


First Moment about 
Plap Hinge (M)B 


Second Moment about 
Shaft Orientation 


eA es Coupling 
Angle (Delta 3) 


Drag Increment 
Lat. Cyclic Pitch 
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DX (15) 

Tw (15) 
XMASS (15) 
Cc 

DPSI 
FTRL 

BIN 

BPIN 
BPPIN 
XITLIM 


RL 

RD 
XLTOL 
XDTOL 
PMOM 


SMOM 
AG, BGL,CG 
TD3L 


DELD 
A1S 


DIMENSION REQUIRED 


FPS YES a 
FT YES a 
FPS YES i 
SLG/CUFT YES 
- YES 

KTS 89, 90 

- YES 

- YES 

DEG 92 

SLG/PT 78,79 

PT YES F 
DEG “se = 
- ** = 
RAD shad 5 
RAD/SEC  *™ 7 
RAD/SEC##24** 

~ me 

LB YES q 
LB 95 : 
LB +4 7 
LB + 4 
SLG-FT 38-52 , 
SLG-SGET 38-52 

DEG ** 

DEG = 

DEG Hm 5 


a 


86 
87 
88 
a9 
90 
91 


(92 


93 


94 
95 
96 
97 


100 


101 
102 
103 
104-105 


106 
107 


108 
109 


110 
141° 
112-113 


414: 
115 
116 


117 
118 
119 
120 


Long. Cyclic Pitch 
Collective Pitch 
Inflow Ratio 
Advance Ratio 

VEL Control 
Iteration Output 
Linear Twist 


No. of Blade 
Segments 


Climb Rate 
Flat Plate Area 
Thrust Option 


Flapping Re-Use 
[Indicator 


MO Iteration 
Tolerance 


Beta Tolerance 
Beta* Tolerance 
Spar Segments 


Second. Harmonic 
Controi Inputs 


solidity 


Symmetric 
Airfoil 


Spring Constant 
Damping Constant 


Lift Only Option 
Shaft Angle 


Desired .A1 and 
B1 Flapping 


Ai, 81 Tolerance 
Tangent Delta 3 


Phase Angle for 
Delta 3 


Induced Velocities 
Induced Velocities 
Not Used 

Input Spar Data 
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B1S 
[75 
LAMBDA 
MUL 
OIN 
SKIPIN 
TWIST 
XNSEG 


RCFPM 
FPAREA 
TOP 
PCNV 


ABIT 


BTOL 
BPTOL 
SPAR 
A2S, B25 


RSL 
SYM 


SFH 
FDUP 


ALOPT 
ALL 


RA1S 
RBIS 


TOLAB 
TD3B 
PHD3B 


LAML 
OVL 


BSPL 


PT-LB/RAD 


FT-L3/ 
RAD/SEC 


% * 


xs 


rT tk, 


121 


122 


123 


124-125 
126 
127-137 
138 


139-141 


142 
143 
144 
145 
146 


147 
156 
157 
158 
159 
160 
161-175 


Azimuthal Printout 
Indicator 


Minimum Lift 
Curve Slope 


Iteration Gain 
Factor 


Not Used 
Pre-Coning Angle 
Not Used 


Bab Moment Inplane 
Aero Forces 


Aircraft Yaw, Roll 
and Pitch Angular 
Velocities 


Cs Staticn 


Rotor Center Staticn 


CG Waterline 
Rotor Waterline 


Aircraft Lateral 
Velocity 


Spar Symmetry 
Airfoil Tables 
Aizrfoil Tables 
Tip Sweep 
Sweep Station 
Rotor Moments 
Rotor Blade 


Airfoil Data 
desi gowente 


FSCG 
FSR 
WLCG 
WLMR 
VELY 


SYMSPR 
HIALPA 
LOALFA 
TIPSWP 
TPSYWST 
TRIM 
RD 


DEG 


RAD 


RAD/SEC 


INCHES 
INCAES 
INCHES 
INCHES 
KT 3 


** program has automatic default value. 


cn 


*% 


*% 
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GRE INPUT DEPAULT VALUES 


ITEM 
NO 


68 
69 
73 
16 
17 
81 
85 
86 
a7 
88 
91 
93 
97 
101 
102 
114 
121 
122 
123 
156 
157 
159 
161-175 


PROGRAM VARIABLE 


DPSI 
FIRL 
XITLIS 
XLTOL 
XDTOL 
BGL 
A1S 
B1s 
T75 
LAMBDA 
SKIPIN 
XNSEG 
PCNYV 
BTOL 
BPTOL 
TOLAB 
PPSI 
ATEST 
IGc 
HIALPA 
LOALFPA 
TPSWST 
RB 


toa epee 7 Are Re CE tee A gee “i 


DEFAULT VALUE 

15.0 

15.0 

15.0 : 
100. | 
50.0 
90.0 


1.0 


* 0. if ABIT, item 100, is non-zero 


** -50 is TOP, item 96, iS non-zero 
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H. CASE INPOT DATA 


Many of the case inputs are self~explanatcry by their 
name listing alone, however, many are not. This section 
will explain the input variables and case options available 
to the user. 

The program enters all the variables into a 200 element 
array called V({I). Prior to entering case data, the prograa 
will automatically do two things. It will initialize the 
array V(I) to a value of zero. It will assign the default 
values listed in the previous section to those particular 
variables. 

The V(LI) array is associated with the variable names by 
equivalent statements. The user needs only to enter values 
for the variables that are different from the default or 
initialized values. In a multiple case computer trun, the 
user need only enter variables that are different from the 
preceeding case. If no new value is enter for a variable, 
the value from the previous case is carried over. 


1. Item 6 -_ Velocity =_VEL 


The computer program will accerft forward velocity in 
one of two ways. The user will input either VzL, item 6, in 
knots or advance ratio MUY, item 89. The value of UIN, item 
90, determines which variable will be used. If UIN is zero, 
VEL will be used. If UIN is non-zero, MUL will te used. If 
VEL is used, the program calculates the advance tatio by the 
following expression. 


UES 
jb = e (v/ ae) x \ 


If MUL is used, the program will calculate the flight 
velocity by the following expression. 


Me 20 es 


} 
NU = ‘a (ge aN YE i 


2. Item _7_-_Hinge Offset -_ER 


The offset ratio of the flaf hinge, E/k, is the 
distance from the center of the rotor shaft tc the vertical 
flapping hinge, normalized by the rotor radius, item 2. 


3. Items 8-22 - Delta X_- DX 


Delta xX is the non-dimensionalized width of each 
individual blade segment starting with segment number one. 
There may be up to 15 segments entered. The number of 
widths entered here must equal the value of Item 93, XNSEG. 
XNSEG is the number of segments intc which the blade is 
divided. This can range from two to fifteen. It is 
recommended that a value of ten or more ke used for XNSEG. 
If XNSEG were equal to 12, values of Delta X would be 
entered for items 8 to 19 and no values would be entered for 
items 20 to 22. The sum of ER plus the summation of the 
Delta X's must equal one. ER is the non-dimensional width 
between the rotor shaft and the flapring hinge cffset. Item 
number 8, which is the first segment width, represents the 
Width between the flapping hinge offset and the point where 
the actual rotor blade airfoil begins. This area is known 
as the spar or cut out segment if no spar data were enterad. 
If item 103, SPAR, is zero, the program will assume that 
this first section creates no lift and has the drag 
characteristics of the first inputed airfoil data section. 
Since this area experiences relatively low dynamic pressure, 
the calculations in this segment do nct have an appreciable 
effect on the outcome of the program. The area between the 
shaft and the hinge offset, £8, produces neither lift nor 
drag in the program's calculations. The tast segment is 
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considered a tip loss factor segment. The lift is assumed 
to be zero, but drag is calculated in the normal manner. In 
previous runs, the width of this section has been set equal 
to three percent of the rotor radius, or Delta X equal to 
03. 


4. Items 23-37 - Twist - Tw 


es ee Se re a Se ee wre SS a ee ee 


The program has two options fcr entering geometric 
twist into the calculations. If the blade has linear twist, 
item 92, TWIST, can be used. If the twist is non-linear, 
the twist can be entered in items 23-37, Td. If a number is 
entered for item 92, the program will assume linear twist 
and ignore ail values entered for TW. The local twist at 
the center of each segment can be entered starting with item 
23. tf the blade contains ten segments, items 23-32 would 
he entered and no values for items 33-37 would be entered. 

A word of caution is necessary regarding the linear 
twist option, item 92. The twist is considered to be zero 
at the 75 percent chord point. The twist is calculated 
assuming that the twist starts at the rotor shaft and varies 
linearly out to the rotor tip. If the actual roter blade 
airfoil started at the 25 percent radius point with a twist 
value of -9 degrees, a value of linear twist equal to -12 
degrees would have to be entered in crder for the correct 
twist distribution to be calculated by the program. 


5. Items 38-52 - XMASS 


ee ee a ee ee ee 


The program provides two methods fcr entering the 
local mass density or the moment of inertia information. 
The individual mass density for each section can be entered 
in items 38-52, or the Pirst and Second Moment, FMOM and 
SMOM, about the Flapping Hinge, cz. be entered in items 78 
and 79. If data are entered for item 79, the Second Moment 
about the Flapping Hinge, the progran will use the 


information provided by items 78 and 79. If variable 79 is 


54. 
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equal to zero, the program will calculate the First and 
Second Moments for items 38-52 and ignore any value entered 
for items 78 and 79. 


6. Items 53-57 -_ Chord 
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The program has two methods for entering local chord 
data. The local chord at each segment can be entered 
starting with segment number one. If the chord is a 
constant chord, the amount of data to be entered can he 
reduced to only items 53 and 54. If items 53 and 54 are 
equal, the program assumes that the chord is constant 
throughout the radius and will ignore any information 
entered in items 55-67. Therefore, for constant chord, only 
enter values for items 53 and 54. If the rotor solidity is 
not inputed in item 106, the program will compute solidity 
as follows. 


DPSI, or Delta PSI, is the incremental azimuthal 
value by which the program advances in its blade flapping 
and force summation coutine. This number must divide evenly 
into 360 degrees. A minimun value of five degrees is 
permitted. It has been found that for most cases decresing 
the value of DPSI below 15 degrees does not improve the 
accuracy but does increase the computational time. As an 


example, the rotor horsepower required for one particular 
run was 1096 RHP for DPSI equal to 15 degrees, 1095 RHP for 
DPS equal to 10 degrees and 1097 RHE for DPSI equal to 5 
degrees. DPSI has a default value of 15 degrees. 
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8. Item 69 ~- FIRL 
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FTRL is the maxinum limit on the number of times 
that the program will enter the flapping iteration routine 
in.search of a steady state flapping solution, The program 
has an automatic default value of 15 iterations. The 
program usually arrives at a steady state flapping solution 
within three to four iterations. The only method by which 
the user can actually determine the number of flapring 
iterations required would be to use one of the debug output 
options. However, these options will create a huge amcunt 
of output and the user is cautioned abcut their use. 


9. Items 70-72 BIN-BPPIN 
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The program has three options on how to assign the 
values for the inital flapping angle, velocity, and 
acceleration at the PSI equal zero position. The user may 
{i) input the values, (2) have the program itself calculate 
initial values, or (3) accept the default values of zero. 
Option three is the most used option. It has been 
discovered that there is little or no difference in solution 
time between option two and three. Variable number 97, 
PCNV, controls which option is to be used, If PCNV is 
non-zero, the program will use either option one or three, 
The program initially assigns values of zero to BIN, BEIN, 
and BPPIN before the initial case data are entered. PCNV is 
assigned the default value of one. Therefore, with no 
values entered by the user for itess 70-72 and 97, the 
program will start with an initial value for flapping angle, 
velocity and acceleration at the PSI equal zero position of 
zero. 

If the wser sets PCNV equal to zero, the following 
formulas are used for determining the inital values. These 
formulas are derived from page 194 of kef. 4, 
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The term ¥ is referred to as the Lock number where 
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In a run where more than one case ig executed at a 
time, a non-zero value for PCNV allcws the previous case 
values for flapping angle, velocity and acceleration tc be 
used as the initial estimate for these variables in the next 
case. If PCNV is equal to zero, the inital values will be 
Calculated by the above fecrmnilas, Most casea ara cun by 
accepting the default valua of one for PCNY. 


10. Item_73 -_XITLIM 


After the program Calculates a steady state flapping 
Solution for its estimated values of inflow ratio, pitch 
angle and shaft tilt angie, the program determines thu 
forces and moments g>nerated by this solution. If the 
forces and, optionally, moments do not meet the reyguired 
amounts entered by the user, the program will cslculaty now 
values to re~enter the flapping routine. KXITLIM determines 
the maximum number of times the program will compare the 
calculated values to the dagired valueg, The prograg has a 
default value of 15 for XITLIM. A majority of the 
soluticns require approximataly five iterations to converge. 
Once the XLTLIM limit 4s exceeded, tha program will stop and 


printout an "Excegdad Limit" statamant. The frograg 
automatically printa as autpul, the nughbes cf Mayor 
Iterations dt uses Sn calculating the solutvten. 


The sign of XITLIM controls the Short Iteration Option. 
If XITLIM is a negative number, the program uses this 
option. A complete describtion can be found in the section 
entitled Main Iteration Opticns. 


The required lift, RL, should equal the actual 
weight of the helicopter that. the rotor system is 
supporting. 


12. Iken_75 - Required Drag 


The required drag can be entered in one of two ways. 
It can be entared directly in item 75 or it can be 
calculated by the program by entering the aircraft's flat 
plate area in item 95, FPAREA. The drag is equal to the 
flight path force that the rotor system must overcome to 
sustain level flight at a certain velocity. Items 75 and 92 
can be entered either as positive or negative numbers and 
the program will provide the corract fcrward flight solution 
by assigning the proper sign value internally. If a valve 
de antered for PPAREA, the program will ignore any value 
aseigned to AD, The drag for FPAREA is calculated as 
follows. 


D=t eV" C FOAREA) 
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If iten 92 is not entered, the user must supply the 
value foc the required drag by using the above fermula. 
There are certain options for which the program igncres or 
does not iterate on drag. An example of this would be when 
the program 1s used to estimate wind-tunnel test results, 


If these options are desired, seo variables 96, TO, and 
110, ALOPT. 


at 


13. [tems 76-77 - Tolerances 


XLTOL and XDTOL are the lift and drag tolerances. 
fhe program will iterate until both the lift and drag are 
within the tolerances given by variables 76 and 77. The 
program haS automatic default values of plus and minus 100 
pounds for lift and 50 pounds for drag. 
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Information regarding the use of FMOM and SMOM can 
be located in paragraph five on XMASS. 


15. Items 80-82 - Shaft Axis 


The gravitational or weight vector must be oriented 
to the shaft axis of the rotor system. Figure 9 shows the 
positions of these angles. The shaft can be oriented in any 
desired direction. The program will automaticaliy assign 
the proper values for normal helicupter flight. The default 
values for AG, BGL, and cCG are 0, 90, andoO degrees, 
respectively. This orients the shaft in the vertical 
direction for normal flight. If any other orientation is 
desired, the user must enter the appropiate values for items 
80 to 8&2. 


16. Item 83 - Delta 3 
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Delta three inputs are controlled by variables 83, 
115, and 116. These are TD3L, TD3B, and PHD3D, 
respectively. If the flapping hinge is connected in such a 
Manner as to cause the blade to change pitch due toa 
flapping, this is referred to as a Delta Three hinge. Item 
83 is the pitch-flap coupling angle; 115 is the tangent of 
the Delta Three Bar; and 116 is the phase angle for the 
Delta Three Bar. The program reduces the pitch angle ata 
particular azimuth and segment by the quantity TD3 where 
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Figure 10 Shaft and Control Axis Diagram 
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[D3 = TD3L ¢ TD3B*sin(PSi + PHD3D) 
Normally, variables 63, 115, and 116 are zero. 


17. IJten 84 - Drag Increment 


The drag increment is a value of delta CD that is 
added to the value of CD obtained from the airfoil input 
data decks. This is added as a roughness factor that 
naturally occurs on blades that are used on production 
aircraft. A value of .002 is normally used. The value is 
not added to the drag calculations for spar data. 


18. Items 85-86 ~- Cyclic Pitch 


The lateral and longitudinal cyclic fitch is 
controlled by four variables. The program allows the user 
to enter both first and second harmonics of cyclic inrut. 
The variables A1S and B1S, 85 and 86, contrel the first 
harmonic inputs, The variables A2S and B2S, 104 and 105, 
control the second harmonic inputs. ‘the A's correspond to 
the lateral inputs and the B's correspond to the 
longitudinal inputs. The program has automatic default 
vales for A1S, B1S, A2S and B2S of ~1.2, 7.53, QO and QO 
degrees, respectively. The user may enter different values 
4f desired. Unless other oftions are indicated, the program 
will keep these cyclic values constant throughout the run 
and will vary inflow ratio and pitch angle to obtain a final 
solution. The program will automatically calculate the 
postion of the shatt axis by the momentum theory and will 
use the value of B1S to deternine the position of the 
control axis. Figure 10 demostrates this fact. 

There is a different option available which allows 
the program to seek specific values of longitudinal and 
lateral flapping. This reguires the use of variables 100 
and 112 - 114. If this option is taken, the values of A15S 
and BIS are sat aqual to zero initially and will be changed 

i 


by the computer in its iteration of the required flapping 


Sl. 


angles. The program also has an cption to remove hub 
rolling and pitching soments. Once a solution is obtained 
that meets the lift and drag tolerances, the values of A1S 
and B1S are varied to reduce the moments. This option is 


controlled by variable 160. Normally, runs are made with 
the user not inputing values for iteas 985, 86, 104, 105, 
100, 112, 113, 114 and 160, 


19. Item 87 Pitch Angle 
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Variable 87 is the initial value of the collective 
pitch at the 75 percent radius station at PSI equal zero 
azimuthal station. The program has a default value of five 
degrees. This value is used only to initiate thc: progran. 
The program, under the normal run option, will vary this 
value in the process of iterating for a convergent solution. 
There is an cption where the pitch angle remains fixed as in 
a wind-tunnel test. This is the TOP oftion, variable 96, 


20. jtem 88 - Inflow Ratio 


The variable LAMBDA controls the inital estimate of 
an uniform inflow. Since the equations of the program are 
done in a gyrocopter mode, inflow is negative when air flows 
down through the rotor. This is the normal ferward flight 
mode. The program has a default value of -.02 for LAMBDA. 
The program will iterate on LAMBDA in its normal iteration 
rountine. The ALOPT, variable 110, oftion will hold LAMBDA 
constant. 


IO SE ED A OD | EP ED iD OE a EP 


Information regarding the use of UCL and UIN can be 
found in paragraph one on Plight Velocity. 
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Information regarding the use of SKIPIN can be found 
in the section of Case Optional Output Indicators, 
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Information regarding twist can be found in 
paragraph four on local twist. 


24, [tem _93.-_XNSEG 


Information regarding XNSEG can be found in 
paragraph three on Delta X. 


25. Item 94 - Rate of Climb 


The program can be made to calculate a complete 
solution for any given rate of climb or descent. Climk or 
descent rate must be entered in units of feet per minute, 
with positive values for climbs and negative values for 
descents. The program assumes a u,iform down or up flow 
across the entire rotor surface equal to the rate of climb 
Or descent. This value is added as an incremental 
correction into the calculations of UP and will effect PHI 
and angle of attack, 


26. Item 95_=_PPAREA 


Information regarding the use of PPAREA, flat plate 
area, can be found in paragraph twelve on Required Drag. 


27, [tem 96 -_ TOP 


The TOP option is one of the wind-tunnel opticrs. 
If TOP is a non-zero number, the program iterates to obtain 
the required lift of variahle 74, but will ignore the 
required drag Of variable 75. Item 87, tha collective pitca 
at the 75 percent radius at PSI equal zero, will te held 
constant. Item 88, the inflow ratio, will be varied in the 
major iteration routine. A non-zero value of TOP will 
result in a value of ~50 for item 122, ATEST. ATEST is the 
minimum accepable value for the lift curve slope when option 
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96 or 110 is executed. If non-zero values for both TOP and 
ALOPT are entered, the program will do the TOP option. 
Shaft angle, item 111, must ke inputed by the user. 
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Information: concerning the use of PCNYV, the Flapping 
Solution Re-Use Indicator, can be found in paragraph nine on 
Initial Plapping Conditions. 


Information concerning the use of PRINT, the 
program's main output indicator, can be found in the section 
entitled Case Optional Output Indicators. 


Item 99 is the End of Case signal card. It is the 
last data variable that will be entered for each case. If 
XEND is a negative number, the program will stop after it 
determines a solution for that particular case. However, 
since an infinite number of cases can be entered for each 
computer run, XEND also tells the program if there are nore 
cases to go. If XEND is equal to 2.0, the program will 
asSume that tne next case will begin by reading in new 
airfoil data. If XEND is any other fositive real number, 
the program assumes that the next cas¢ will use the present 
airfoil data and will enter only case input data and any of 
the options which normally follow the case input data. Each 
time that the variable XEND is entered, be especially 
careful to follow the format for NNUM for this card. NNUM 4 
ds the pumber of inputs per data card. NNUM must be a 
negative number for this XEND card. It is this negative 
3ign on NNUM which actually keys the computer to stop 
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Information regarding the use of ABIT can be found 
in paragraph 40. 


32. [tems 101-102 - BToL 


In the blade flapping routine, the program searches 
for a steady state flapping solution for the given 
conditions of inflow ratio, pitch angle, and cyclic input. 
The program compares values of flapping angle and velocity 
at the PSI equal zero azimuthal positicn on each revolution. 
If at this position, the difference between the n-th and the 
(n - 1)th revolution values for flapping angle and velocity 
is less than BTOL and BPTOL, respectively, the program 
assumes that it has determined a steady state solution. 
BTOL and BPTOL have default values of 0.000001 radians and 
radians per second, respectively. If other values are 
desired, the user may enter these values for items 101 and 
102. 
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The number of blade segments using spar data can be 
inputed in item 103. If ne spar data are available, no 
value for SPAR should be entered. For this case, the 
program assumes that the first segment is the area between 
the flapping hinge and the point where the actual airfoil 
begins on the rotor blade. This area is also referred to as 
the "cut out" seqment. In this case, the program assumes 
that this cut out area produces zero lift and uses CD 
information from the first aixfoil section for drag 
calculations. If a non-zero number is entered for SPAR, 
spar airfoil data must be availabie. Case input variable 
120, BSPL, controls the spar input option. If spar data are 
to be inputed, BSPL should be assigned a non-zero value. In 
a multiple case run, where spar data are initially entered, 
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the variable BSPL is set equal to zerec as soon as the spar 
data are entered. If the program did not do this, the user 
would have to enter a zero for BSPL for the next case if no 
new spar data are to be entered. If in a multiple case trun, 
the user decides to enter new spar data, the variable BSPL 
must be assigned a non-zero value for that particular case. 


34. Items 


Information regarding the use of the second harmonic 
control inputs can be found in paragraph 18 on Lateral and 
Longitudinal Cyclic Inputs. 


Information regarding the use of RSL, retor 
solidity, can be found in paragraph six on local chord. 


36. Item 107 _- SYM 


SYM is the non-symmetrical airfoil input control. 
If the user assigns a non-zero value for SYM, the program 
will assume that all blade airfoil data are non-symmetrical. 
The user must enter values for CL and CD for the complete 
range of angles of attack from -180 to +180 degrees. If the 
value of SYM is zero, only tabular values from zero to +180 
degrees need to be entered. The above holds true also for 
variable 147, SYMSPR. SYMSPR applies to the spar data 
exactly in the same manner as SYM applies to the airfoil 
data. 


Values for the spring constant, SFH, and damping 
constant, FDMP, about the flapping hinge can be entered if 
known. These variables can be entered to simulate a 
hingeless rotor system or a system with flapping springs. 
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38. Item 110 - ALOPT 


This is one of the wind-tunnel options. If ALOPT is 
a non-zero input, the program will iterate to obtain the 
lift required of variable 74 but will ignore the required 
drag, variable 75. In this option, variable 87, collective 
pitch, will be varied, but not variable 88, inflow ratio, in 
the program calculations. This is the opposite of the TOP, 
varialbe 96, option. 

If a lift curve slope less than ATEST, item 122, is 
calculated while using this option, the program will stop 
and produce the following message. "Stall Criterion has 
been violated -- will go to next case, if any." Item 111, 
the shaft angle, must be inputed., Jf non-zero values are 
entered for both TOP and ALOPT, the program will do the TOP 
Option. 
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Item 111, the shaft angle, must be inputed whenever 
the TOP or ALOPT options are used. 


40. Items 112-114 = RAIS 


Yf variable 100, ABIT, is non-zero, the program will 
iterate the blade flapping solution in an attempt to obtain 
the desired lateral and longitudinal flapping angles 
indicated by variables RA1S and RB1S, respectively. The 
program will iterate to the accuracy indicated by TOLAB, 
item 114. Item 112 is RA1S and item 113 is RB1S. 


Items 115 and 116 are the tangent and phase angle of 
a Delta 3 Bar. [Information regarding the use of TD3B and 
PHD3B can he found in paragraph 16 on the Pitch-Flap 
Coupling Angle. 
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The program allows the user tc induce a velocity of 
any form onto the rotor systen. This is done in a harmcnic 
series of the form (A0 + A1*cos(PSI) + Bt*sin(PSI) + 
A2*cos (2*PSI) + B2*sin(2*PST) + ....) for each segment that 
the blade is divided into. If variable 117, LAML, is a 
non-zero number, the program will enter the harmonics of the 
induced velocities. During the Read routine, LANL, will be 
assigned a value of zero. Therefore, for each ‘case where 
different values for the harmonics are desired, LANL will 
have to be set toa non-zero number. Variable 118, UVL, 
controls the use of the induced velocities. If UVL, is 
zero, the induced velocities will all te set equal to zero. 

A short example will now be given on the use of the 
control variables ina multiple case run. Assume that no 
harmonic induced velocites are desired for the first case. 
The user would make no inputs for LAML and UVL. For the 
second case, assume that harmonic induced velocities are 
desired. LAML and UVL would be set equal to a non-zero 
number, The harmonic variables would follow the case input 
data for this particular case. For the third case, it is 
desired that the same induced velocites be used. The user 
woul? not enter any values for LAML and UVL since (1) LAML 
has been automatically set to zero and hence ns new harmonic 
data will be entered and (2) UVL is still equal to a 
non-Zery nugder. It is desired for case four to use no 
induced velocities, The user now would enter the value of 
zero for UVL and program will zero out all the induced 
velocity variables. ‘ 

This paragraph will describe how to use the variable 
induced velocity option. The inputs Ail, Bi, A2, B2, «4. 
are numbers in units of feet per second. If the direction 
of the velocity is down, a negative sign is associated with 
the A's and B's. A negative sign will decrease the angle of 
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attack for an element originally at a positive angle of 

- attack. It will increase the amount of negative angle of 
attack for a blade element originally at a negative angle of 
attack. The values of A‘s and B's do not effect the 
uniform inflow ratio, LAMBDA, that the program iterates 
upon. The format for using the variable inflow velocity is 
as follows. 

1. A value for NHARM is entered in I3 format. 
NHARM is the maximum degree of the harmonics 
that is to be used. If the maximum degree used 
is A3*cos(3*PSI), then NHARM is three. 

The next three paragraphs are repeated for each 
blade segment. 

2. The value for AO is entered in £E15.6 format. 

The IBM 360-67 at the NPS will accept F15.0 
format. 

3. The values for Al, A2, A3, «es, Up tO NHARM are 
entered on the next data cards in format 
5E14.6, 

4. The values for Bi, B2, B3, «ee, Up to NHARM are 
entered on the next data cards in format 
S5E14.6. 

5. Three cards are reguired for each blade segment 
whether or not the values are equal to zero. 

If the user has 15 blade segmen’ , a minimum of 
45 data cards are required. 
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Information concerning the use of BSPL can be found 
in paragraph 33 entitled Number of Spar Airfoil Data 
Segments. 


PPSI represents a delta PSI for printout purposes. 
PPSI has a default value equal to DPSI, variable 68. PPSI 
Shall never be a smaller increment than the incremental DPSI 
used to calculate the soluticn. 
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tem_122_-_ATEST 


ATEST is the minimum acceptable value of the lift 
curve slope when option TOP or ALOPT is used. Tf no value 
is assigned, LTEST has a default value of 5 (1/rad). When 
the TOP option is used it has a automatic value of <-50 
(1/rad). The lift curve refers to the increase in rotor 
lift with the tilting back of the tip path plane. A check 
on it for a minimum is for ccnvergence purpose only. 


Icc is the Iteration Gain Control factor. If the 
major iteration fails to converge, choosing a fractional 
value for IGC can greatly speed convergence. This may he 
especially helpful when parts of the rotor are in stall. 
The amounts Ly which the convergence algorithm changes the 
independent variables is muitiplied ty IGc. Setting item 
91, SKIPIN, equal to zero may help the user decide if this 
IGC option might be useful. 


47. Item 126 ~ PCR 
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48. Item 138 - INPL 


Input 1.0 for INPL in order to remove huh monent 
inplane aerodynamic forces from the calculation of 
aerodynamic pitch and rull moments about the shaft axis. 
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49, [tems 139-J4l_- psis 
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Variables PSYS, VHIS and THES rapresent aireraft 
angular velocites. The ajrceraft can be givan any argular 
valocity in yaw (BSIS), pitch (FHTS), and roll (THPS) in 
radians per second by the use of theso variables. They 
effect the calculations of UP, UT and UL. fFaragraph $0 


contains additional information. 


59. Lheme 142-199 -c¢ 


If the angular velocities, ft¢eny 139 ~ 141 are 
entered, they are used in the calculaticns of UP, UT and UR. 
Tf no values for items 142 = 145 are entered, the program 
AUSUVOS that the retor wybtem rotates due to the angular 
yalocities about the center of the a#haft. If values are 
entered for ditegn We - 145, ths calculations will easaspuge 
that the cntite roto, shaft de cuotatinug about the center of 
gravity. Varsables 142 - 949 ara ocG, PSMK, Wicd and Wher, 
¥SCG de the longitudinal Cy position whiso FLUMP de tho 
Longitudine) position of the malt rotor, #LG4 iB the GG 


‘Waterling stotfion while WLAP de the main toterta watsrline, 


OL Verticad position, All of theae values wast be cutelad 
fu units of inches, 
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O2- [sna 142.2 Sxbeke 


If the spar date are syomatrical, 40 not enter a 
value for SYMSPR,. The proygvag wil) only use valwon for CL 
and Cb between 4are to 4189 deyices, Ff the mpat data area 
non-symmetrical, enter any non-zeto valn@ for SYMobih. The 
program will veqguite Yalugs from -180 to #180 deg: ave. 
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53. liesg 126-127_-_HIALFE 


In order for the program to properly calculate the 
region in and around the reverse flcw region, the valu: s of 
CL and CD are required to be Known at high and low angles of 
attacks approaching 140 degrees. However, in order to save 
the user from naving to enter a whola range of angles of 
attack for all Mach numbers, the program has an option where 
for angles above HIALFA, 156, and below LOALPA, 157, the 
Mach number is set equal to zero for table lookup purposes. 
The user ds only requircd to enter large angles for the 
first two mach nunber tables. HIALFA and LOALFA have 


54. [tems 128-159 _-_Tip Sweep 


The GRP prcqram will properly calculate the airflow 
sYGep angle, Ui, UK, and pitch angle of a swept tip airfoil. 
TIPSWP js the amount of SWuep of the tip measured in 
degrees. TeSastT ts the blade seqment number at which the 
gveep begins. The program assumes that the remaining 
outboald sagments starting with TPSWST are swept the number 
of degrees Lndicated by TIPSWP., UT and UR are modified for 
thes@ segmants as shown is Figure 11. 


99- [hqp.109_-_TEIA 


If TRIM dis a non-zero value, the frogram will 
attempt to adjust Als and B1S in order to reduce the rolling 
and pitchiny woents to less than plus or minus 190 pounds. 
ye owhll tduet obtain a solution which will satisfy the 
tegudred lift, and diag, then it will adjust A1S and B18 to 
reduce the moments. 
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Variables 167 - 175 assign the proper airfoil data 


to blade segments one through fifteen. The program will 
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accept up to five airfoil data decks. The RB array is 
initialized with the value of one for all 15 segments. If 
only one blade airfoil deck is used, there is no need to 
enter any values for RB. Only enter values for the segments 
which will use airfoil data sets two through five. If no 
Spar data are entered, SPAR equals zero, Segment number one 
is considered the cut out segment and is assigned the value 
of RB(1) equals one for calculation purposes and the value 
of RB(1) equals zero for printout indentification. In this 
case, segment one produces only drag but no lift. If SPAR 
is greater than zero, RB(I), I = 1, SEAR, will be assigned 
the value of zero for printout indentification. These 
segments will use the spar data entered and will calculate 
both lift and drag. 


I. CASE INPUT FORMAT 


Below is located the format that ic used to input cuse 
data to the program. All input data cards are cf the format 
(12, I4, 5F12.0). The input cards contain seven fields 
Which are called NNOUM, NLOC, C(1), €C(2), €(3), C(4) apd 
C(5). An example of a data input card is as follow, 

5 1 700. 26.8 1148.6 002246 4.0 

1. NNUM ds the number of inputs on this card, 
where C(N) are the inputs. NNUM must appear in 
either cclumgn one or two. NNUM has a oininum 
of one and a maximum of five. In the ahova 
example NNUM is five, 

2. NLOC is the item or variable number of input 

C(1). This refers to the item nuwhers that ara 
found in the section on Casa Input Listirgnu. 
NLOC is in 14 format and must be cight 
justified in columns four through oix. In the 
example NLOC refers to item number one, the 
Rotor Tip Speed, 
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3. (1) through C(5) are the values corresponding 
to the input items NLOC through NLOC + NNUM. 
Each value must contain a decimal point and be 
dn columns 7 = 18 for €(1), 19 - 30 for C(2), 
371 = 42 for €(3), 43 - 54 for C(4) and 55 - 66 
for c(5). In this example, values are entered 
for variables one through five, the Tip Speed, 
Radius, Speed of Sound, Air Density and Number 
of Blades. 

4, Omision of NNUM will cause program termination 
with an error explanation statement. NNUM 
greater than five will cause unknown problems. 
Omission of NLOC will cause the present card's 
values of C(N) to he entered into the items 
indicated by the previous card's NLOC. Failure 
to right justify NLOC, or NLOCC greater than 
200, will cause unknown problens. Failure to 
properly locate correctly any input value 
within its own field on the card will cause 
exrors in both that jnput and tne number whose 
field it encroaches on. 


J. CASE OPYIONAL OUTPUT INDICATORS 


Tha program has two variables that control the output 
printout, variable 91, SKIPIN, and 98, PRINT. The prograo 
will automatically produce a printout of the case input data 
for each case and a one=~page summary of the inital 
conditions and iteration limitations thea user placed upon 
tha progcam. It will also produce a one-page summary of the 
raaulting forces, goments and calculated rotor horsepower 


tor the ginal converged solution. The user can also receive 
an echo printout of the airfoil and spar data decks. If 
this acho printout ds desired, the user will find the 
correct printout indicators described in sections C and D 


entitled Blade and Spar Data Deck Requirements. 

The main optional printout variable is variable 98, 
PRINT. PRINT can be a number from 1 to 1,111,111 depending 
upon the option desired. If no value is enter for PRINT, 
the user will receive the printout described atove. Below, 
is listed the PRINT Options. If, for example, PRIKT is 
assigned a value of 111, the user will receive printout 
options 1, 10 and 100. 


OPTIONAL OUTPUT INDICATORS 


Option 1 Angle of attack, Mach number, 
section lift and drag 
coefficients, inflow angle, lift, 
and sweep angle at each azimuthal 
position for each radial blade 


segment. Only for converged 
flapping solution. 

Option 10 Converged flapping angle, rate, 
and acceleration at each 


azimuthal positicn. 

Option 100 Converged integrated forces on 
blade at each azimuthal station. 

Option 1000 Harmonic analysis of blade forces 
for converged case. 

Option 10000 Harmonic analysis of air loads 
for ccnverged case. 

DEBUGGING OR TRANSIENT OPTIONS ; 

Option 100000 Transient flapping angle, rate, 
and acceleration at each 
azimuthal station. 

Opti > 1000000 Option 100009 plus blade 
velocities, angles, Mach number, 
section coefficients, and lift 
for each blade segment at each 
azimuthal station, 
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The user is cautioned that the debugging options 
can give a huge awount of output data. 


The second printout option is the variable SKIPIN, 
number 91. This variakle controls the suamary ferce, noment 
and horsepower output discussed in the first paragraph. If 
SKIPIN is greater than zero, this summary will be printed 
only for the final converged solution, It SKIPIN is equal 
to zero or a negative number, this summary will be printed 
for each loop through the major iteration (force summation) 


routine. SKIPIN has a default value of one. If the progran 


is not converging to a solution, the user can se9 
immediately, with very little extra printout, exactly what 
intermediate solutions the progran is producing by setting 
SKIPIN equal to zero. This may help the usar in deciding 
whether or not to use the Iteration Gain Pactor, IGC, 
variable 123, 


Ke. IBM 360 EXECUTION CONTROL CARDS 


This section idilustrates the control cards required to 
execute the GRP programa using the IBM 360 at the Naval 
Postgraduate School, The program may ba run under OS or 
CP/OMS, There are two ways of running the program undar 0%. 
The first 4s to run the entire program and data through the 
computer at the same time. The second way ie to store the 
main program on a disk as a library program and enter only 
the data through the card reader for each desired case. The 
second method has the two advantages of (1) not reguiring 
the user to enter the entire 1100 card aadn program through 
the card reader for aach run and (2) the amount of CPL time 
required can be reduced since tho main program does not have 
to be recompiled for each cun. the plogiam raquires 
approximat2ly one minute and forty seconds of CPU time to 
compile, The normal run time for aach casa iu approxiuately 
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20 CPU seconds. This can vary with the amount of printout 
data requested. If the program is compiled on the CP/CMS, 
the user will have to request 344K bytes core size on the 
login message. The standard 256K bytes core size is not 
large enough for compiling. However, once the program has 
been compiled, it can be executed within the 256K normally 
available. 

The following cards are required to execute the entire 
program through the card reader at one time. 

Standard Job Card 

// EXEC FORTCLG, REGION. FORT=150K 

// REGION.GO=180K 

//PORT.SYSIN DD * 

Main GRP Prograa 

/* 

//%0.SY¥SIN DD * 

Case Input Data 

/* 

The following two programs are used to reserve space and 
load the program onto a disk. 

Standard Job Card 

// BXEC FYGM=IEPBRIY 

//LOAD DD DSN#S139S HELO, UNIT#3330 ,VCLESER=DISKO1 

// DISE™ (NEW, KEEP) ,LABELSRETPD= 150 ,SPACE=(CYL, (1,1,1)) 

/* 

Standard Job Card 

// EXEC FORTCL, REGION. PORT#180K 

“S/PORT.SYSIN DU * 

Main GRP Prograng 

/* 

// INK. SYSLNOD DD UNIT#3330,VOL"=SERSDISKOI, 

// OSN=@31395. HELO (GRP) ,DISP#SHA 

The $1395 used ahove and helow must ba change to Ss for 
go faculty with the appropiate user number 
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student or F f 
395. The follvewiig catds Gust be ussd to 
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execute the program once stored on a disk. 
Standard Job Card 
//GO EXEC PGM=GRP, REGION=180K 
7/STEPLIB DD UNIT=3330, VOL=SER=DISK01, DISP=SHR, 
// DSN=S1395.HELO 


//¥T06F001 DD SYSOUT=A, DCB= (RECFM=FBA, LRECL=133, BLKSIZE=3325 


//¥PTOSFO01 DD * 
Case Input Data 
/* 


J. Sample Program Output 


This section describes in detail the output available 
from the GRP program. In addition, a sample computer output 
is included with each describtion. The program will print 
up to ten different tables. Seven of these tables are 
optional and are not automatically printed. The ten tables 
are as follows. 

1. Echo Printout of Rotor Blade and Spar Data 

2. Case Input Data Card Listings 

3. Summary of Input Data 

4. Debugging or Transient Information (Options 

100,000 and 1,000,000) 

5. Summary of Forces, Moments and Horsepower 

6. Converged Flapping Sclution (Oftion 10) 

7. Converged Integrated Forces (Oftion 100) 

8. Harmonic Analysis of Z Force (Cption 1000) 

9. Converged Blade Analysis (opticn 1) 

10. Harmonic Analysis of Air Loads (Option 10,000) 

The input variable PRINT, item 98, controls the output 
of items four and six through ten. Items two, three and 
five are always outputed. Item one is contrelled by the 
first data card on the rotor blade and spar section input 
decks. 
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1. Echo Rotor Blade Printout 


The next page contains a partial sample Echo 
Printout of Rotor Blade Input Data. This printout 
illustrates that (1) the printout was requested (WBLADE = 
10.), (2) there are two airfoil decks to be entered, (3) 
there are nine Mach numbers for which CL's are to be read in 
and (4) the remaining portion of the printout is the values 
of Mach numbers, angles of attack and lift coefficients for 


the first airfoil deck. 
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The next Page contains a Sauple Case Input Data Card 
Listing. This is one of the automatic printouts, y-* 


is an 
echo printout of the input data cacds 
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3. Summary of Input Data 


The next page contains the automatic Summary of 
Input Data printout. The following information can be seen 
on this sample output. 

a. The blade was divided into 15 segments starting 
from the hinge offset and proceeding outward. 

b. NO values were entered for the local sass 
desnsity, input variables 38 - 52. Instead, values for the 
Pirst Moment, M(M) = 85, and the Second Moment, SOK ~< 
INERTIA = 1450, about the Plapping Hinge, input items 78 and 
79, were entered. 

c. The X row indicates the calculated centers of 
each segment cxpressed in a percentage of the distance out 
the rotor blade. 

a. The ELADR DECK row indicates that the first 
blade segment was considered 4 Spar segment. Segmanta two 
through five and twelve through fifteen belong to airfoil 
data deck number one, while segments six through eleven 
belong to airfoil data deck nuaber two, 

@. The rest of the information, with one excaption, 
is a summary of the cave input data. The exception im the 
tera THRUST FACTOR, This is the value veed to 
hondimensionalize all the calculated forces in the progran,. 
The TRRUST PACTOR equals @ re CouKy* . Howents ara 
nondiuensionalized by the THRUST FACTOR times the radius. 

f. The progias checks to see 4¢ all the blade 
bey “s, delta X's, plus the distance from the shaft to the 
hinge Ceut, EB/R, 244 up to one. If, on the printout, 
SUM/DX, ¢ P/R Goes not egual one, the user hae maie a 
mistake somauvhere with the delta X's o: in the M/R nugber, 
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4. Debugging of Transient Information 


The next page illustrates the Debugging or Transient 
Printout. This is PRINT Option 1,900,000. The information 
available includes flapping angles, rates and accelerations 
at each az.mutial position with a radial position display of 
UP, UT, U, PHI, Angle-of-Attack, Mach Number, CL, CD and 
Lift Per Inch produced. 

This option is generally outputed only when the user 
is experiencing unknown difficulties with the GRP program. 
The program will output all of the above information for 
every revolution and iteration until a ccnverged solution is 
obtained or the program runs out of allowable computer tine. 
If desired, the variable SKIPIN, item $1, will provide a 
printout of the Porces, Moments and Horsepower Summary after 
each major iteration. 
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5. Forces summary 


The next printout illustrates the Summary of Forces, 
Moments and Horsepower. This page is automatically outputed 
foc the converged solution. A printout of this summary can 
be obtained for each loop through the major iteration 
routine by the use of the input variable SKIPIN, itea 91. 
The following information can be observed from thic saaple 
printout. 

a. The cyclic lateral and longitudinal inputs, AS, 
B1S, A2S an&é B2S are printed at the tcp of the page. This 
exatple indicates that there were no sacond harmonic inguts 
for 425 and £43, which is normal, 

). THETA .75 148 the pitch angle of the blade at the 
76 percent radius station at the PSI equals zero sziauthal 
position, In wost options the program wil) iterate upon the 
values for Thzr2 .75 in dt# search for a convergad solution, 

c. LABE!A refers to the converged value for the 
unifora inflow ratio. If arate of climb of dercent was 
used in the ces, that rate divided by the tip speed wouid 
have to be subtracted f19m or added to this value of LAMBDA, 
respectively. 

G. MU(Y)S and MU(Y)S are the advance ratios in the 
longitudinal and latezal ditections with respect to the 
Shoft axia. 

e. CT, CQ, CH, CL and Cb are the calculated overall 
coefficients of Thrust, Turgue, it Forces, Lift and Diag. Add 
of those ‘items are nondisonsionalized by the value of the 
Thrust Factor. 

f. Lift and Dray forces are calculated with taepect 
ty the tLelative wind axis, Theust and i forces are 
calculated with respect to the control axig. % forces area 
calcuated with tLreppect to the shaft axin, X and ¥ fotcen 
are calculated perpendicular to the shaft axia. 

g. The Equivalent Diag ds the total diay force 
created by the fuselage, flat plate area timuu dynapic 


pressure, plus the profile drag created by the turning rotor 
syYntea. ° 

h. The Equivalent P. A., Or Equivalent Flat Plate 
area is obtained by dividing the Equivalent Drag by the 
dynanic prossure. 

i. Alpha(S) is the shaft axis orientation while 
Alpha(c) 13 the control axis orientation. The program uses 
momentum theory to calculate Alpha(S). Alpha (C) is 
Getermined from the following relationship, Alfha(sS) = 
Alpha(c) + p13. 

J. LAT. OFS. and LONG. DIS, are the lift vector 
offset as a percentage of the rotor radius frem the rector 
@haft to create the prograrg's rolling and pitching goments. 

K. PM and &M aro the calculated aerodynamic pitch 
and coll moments about the shaft axis. The SHEARS Hub Pitch 
and Roll Moments are calculated from summing the 
aerodynamic, inertia and elastic hub restraint moments. 
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6. Flapping soluticn 


The next page contains a sagple printout of the 
converged flapping solution. The first item tc appear are 
the flapping values at the PSI equal zero azimuthal position 
for each revolution prior to the converged ravolution on the 
final iteration. The values for flapping angle and its 
first two derivatives are expressed in radians. 

a The second part of this printout is the actual 
flapping values for the converged revoltion. The difference 
between the flapping angles and rates between the PSI equal 
zero and 360 degree position must be less than the 
tolerances entered in variables 101 and 102. The nusbera 
here are also in radians. The last item is 4 Yourier 
coefficient series for the flapping angle and it ds 
calculated in degrees, All calculations are done in 
respects to the shaft axis, This is PRINT Option 10. 
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7. Force Integration 


The next table is a sample Force Integration Outrut. 
This is PXINT Option 100. The forces are a printout for one 
blade only at a particular azimuthal fosition. CQ, CQL and 
CQD are the Coefficients of Torque, Torgue due to Lift and 
Torque due to Drag. CQ is calculated from CQ = CQD = CQL. 
CX, CY and CZ are all related to forces in the shaft axis 
reference systeng. CMHS is the Coefficient of Pitching 
Mo2ent due to aerodynamic, inertia and hub elastic restraint 
moments about the Shaft Axis and CLHS is the Coefficient of 
Rolling Moment due to these same forces. In the printout 
the (B) character is the number of tlades, which in this 
example is four. SIG is the solidity of the rotor systen. 
MAX B*COD/SIGMHA is a blade stall indicatcr. It is 
calculated by determining the first azimuthal station 
between the PSI equal 180 degree and 360 degree position 
that has a value of CQD greater than the value of CQD at the 
PSI eygual 180 degrees azimuthal position. 
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e. 
4 
: 8. Harmonic Analysis of Z Force 
The following is a harmonic analysis of the forces = 
in the Z or shaft axis directicn. Z force is positive in 
the downward direction. This is PRINT Option 1000. 
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The next four pages contain the sample overall 
summary of events occuring on the blade at each azimuthal 
and radial position. This is PRINT Option 1. Variable 121, 
PPSI, controls the azimuthal intervals that are printed out. 
The output items are as follows. 

a. X - Center location of the blade segnent 

b. ALPHA - Angle-of-Attack 

c. MACH - Local Mach number 

d. CL - Local Coefficient of Lift 

e. CD - Local Coefficient of Drag 

£. PHI - Local Inflow Angle 

g. L(LB/IN) - Lift produced per inch on segment 


2a Ht etal adit abodes Aah 


eet ac at beaded sna ha 


h. Sweep - Sweep Angle of airflow 
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10. Harmonic of Air Loads 
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The next printout contains the harmonic analysis of 
the lift generated per inch on each rotor blade segment. 
This is PRINT Option 10000. The airloads are computed in 
two harmonic forms, both containing terms out to and 
including the tenth harmonic. The forms are 

Lift per Inch = AQ - Ailcosy - BisinY- A2cos27 - 
B2sin2¥ - A3cos3¥ - B3sin3Y ..... = Aldcosl10¥ - 
Bi0sini0¥ and 

Lift per Inch = AQ - Cilsin(w+td) - C2sin(2¥+@) 
~C3sin (3¥+4) ..... - C10sin (10N+4). 

The Ratio column is determined from the coefficients 
in Column C. From this ratio, one can immediately determine 
which airload harmonic is dominant and the relative 
relationship of this to the other harmonic coefficients. 
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IV. GRP SAMPLE ANALYSIS 

The GRP program was executed using data representing a 
relatively new rotor blade. The results were ccmpared with 
results predicted by the blade's manufacturer. The rotor 
blade was an unsymmetrical blade. The rotor radius was 
divided into three sections. Sections one and three vere 
made of the same type airfoil. The blade included a sweep 
tip design. The blade and helicopter configuration analyzed 
are typical of a helicopter that could be used by the J.S. 
Navy in a LAMPS type mission. 

For the analysis, the blade was divided into 15 
segments. The program used the manufacturer's values for 
the First and Second Moment of Inertia about the Flapping 
Hinge, vice local blade mass densities. The GRP assumed 
uniform inflow for all flight velocities . It used a rigid 
blade analysis while the actual blade does have live twist. 
The program was run at five different flight weights, 
ranging from 16,359 to 20,829 pounds. A flat plate area of 
35.8 square feet was used at all speeds. The GRP was run 
for forward flight speeds of 40 to 160 knots at ten knot 
increments. The program was executed at sea level, tropical 
day condition. The manufacturer's predicted rotor 
horsepower was obtained from his Shaft Horsepower versus 
True Airspeed curves and was corrected to rotor horsepower 
by using the manufacturer's Mechanical Efficiency curves. 

The results of the analysis are shown in the next 
several tables. fable VI illustrates a comparison of the 
GRP required rotor horsepower divided by the manufacturer's 
required rotor horsepower, The GRP agreed within an average 
of two percent on the entire range frem 40 tc 160 knets. 
The GRP agreed within an average of one percent for the 
cruise range between 70 and 140 knots. It can be seen that 
between 40 to 60 knots there is amuch larger difference 
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between the two required horsepowers. It is felt that the 
inflow in this region is not uniform as assuaed, but highly 
mixed and irregular. Also, the GRP results were less than 
the wanufacturer's horsepower in the 150 to 160 knot range. 
This area represents the region of top speed for the 
helicopter, and much of the retreating blade is in the stall 
region. Also, it is expected that there is a change in 
fuselage attitude at this high speed, which would increase 
the flat plate area above what was used in the program. The 
GHP program's maximum endurance velocities agreed exactly 
with those predicted by the manufacturer. 


TABLE I WEIGHT = 16359 LBS 


VELOCITY GRP RHP MANUFACTURER'S RHP RATIO (GRP/HAN) 
40 1105 1152 96 
50 1006 1041 «97 
69 960 973 99 
70 955 958 1.00 
80 984 971 71.01 
90 1047 1034 71.01 

10u 1142 1103 1.03 
110 1273 1237 1.03 
120 1438 1396 1.03 
130 1643 1589 1.03 
140 1893 1852 1.02 
150 2300 2216 1.04 
160 2567 2745 ~94 
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| 
TABLE II WEIGHT = 17321 LBS 


VELOCITY. GRP RHP MANUFACTURER'S RHP RATIO (GRP/MAN}) 
40 1185 1247 95 
50 1075 1122 96 
60 1018 1038 98 
; 70 1008 4009 1.00 
E 80 1028 1028 1.00 
90 1086 1078 1.01 
100 1176 1153 1.02 
110 1303 1264 1.03 
120 1465 1445 1.01 
| 130 1669 1633 1.02 
140 1921 1902 1.01 
150 2227 2257 .99 
160 2598 2791 . .93 


TABLE III WEIGHT = 19246 LBS 


VELOCITY GRP RHP MANUFACTURER'S RHP RATIO(GRP/MAN) f 
40 1366 1449 £94 H 
50 1229 1301 £94 
60 1148 1186 97 } 
70 1118 1142 98 ' 
80 1126 1153 .98 ; 
90 1172 1199 .98 

100 1255 1273 -99 
110 1375 1387 .99 
120 1530 1557 . .98 
130 1731 1746 .99 
140 1983 2038 97 
150 2290 - 2440 94 


160 2665 - - ‘ 
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TABLE IV WEIGHT = 19658 LBS 


VELOCITY GRP RHP MANUFACTURER'S RHP RATIO (GRP/MAR) 
40 1407 1496 94 
50 1264 1340 94 
60 1178 1231 -96 
70 1142 1173 97 
80 1149 1179 -98 
90 1193 1221 98 
100 1274 1299 98 
110 1392 Wy 98 
129 1546 1566 99 
130 1747 1773 99 

140 1999 2048 98 
150 2305 2477 ~93 
160 26 82 = - 


TABLE VY WEIGHT = 20829 LBS 


VELOCITY GRP RHP MANUFACTURER'S RHP RATIO(GRP/MAN) 
49 1532 1627 -94 
50 1370 1457 294 
60 1269 1338 .95 
70 1222 1267 .96 
80 1219 1258 .97 
90 1256 1297 297 

100 1330 1375 £97 
110 1445 1477 .98 
120 1598 1638 .98 
130 1798 1866 96 
140 2045 2157 95 
150 2354 2589 91 
160 2747 - - 
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VELOCITY 
40 
50 
60 
70 
80 
90 

100 
110 
120 
130 
140 
150 
160 


AVERAGE 


AV ERAGE 


TABLE VI RATIO COMPARISON 


WEIGHT 
16359 17321 19246 19658 
96 95 «94 294 
97 96 94 94 
99 98 97 ~96 
1.00 1.90 98 97 
1.01 1.00 98 97 
1.01 1.01 98 -98 
1.03 1.02 -99 98 
1.03 1.03 99 99 
1.03 1.01 98 99 
1.03 1.02 99 99 
1.02 1.01 97 98 
1.04 -99 294 293 
94 93 = = 
AVERAGES FOR ENTIRE SPEED RANGE 

1,00 99 97 97 


AVERAGES FOR CRUISE RANGE 
70 - 140 KNOTS 
1.02 1.01 98 98 


1i2. 


20829 


~94 
94 
«95 
96 
37 
97 
Pie De 
98 
98 


- 96 


97 


RATIO 
AVERAGE 


295 
295 
97 
-98 
99 
99 
1.00 
1.00 
1.00 
1.00 
99 
- 96 
«94 


-98 


99 


V. CONCLUSIONS 


The logic and theory used in the GRP program was 
investigated and found to be sound. However, there were three 
discrepancies in the Navy's version of the program that did 
The calculations in the reverse flow section 
the calculation of the chord at the 


require attention. 
on the rotor were incorrect, 
75 percent radius position was incorrect and the original Trim 
Option for reducing moments would not work. All of the above 
discrepancies were corrected. 

Tnree desirable features were added to the GRP program. 
First, the ability to analyze a rotor blade composed of more than 
one airfoil type was added. The program will now accept up to 
five different airfoil data input decks for use in analyzing a 
rotor system. Secondly, the program would only calculate performance 
in level flight. The ability to calculate performance in climbs 
and descents has been added. Lastly, the program will now calculate 
the aerodynamics for a swept tip rotor blade design. 

The results of the sample analysis described in Section IV 
indicates that the program does produce highly accurate performance 
predictions. The averaged GRP rotor horsepower was within two 
percent of the manufacturer's data. The results were within one 
percent when compared in the area of normal cruise flight. The 
GRP, in this analysis, assumed uniform inflow, constant flat plate 
area and a rigid rotor blade. It is felt that while complicated, 
computer-time-comsuming procedures can be taken to reduce these 
assumptions, they are not warranted if the GRP is to be used 


strictly as a helicopter performance prediction program, 
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